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Nonextensive relativistic nuclear and subnuclear equation of state 
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Following the basic prescriptions of the Tsallis' nonextensive relativistic thermodynamics, we investigate the 
relevance of nonextensive statistical effects on the relativistic nuclear and subnuclear equation of state. In this 
framework, we study the first order phase transition from hadronic to quark-gluon plasma phase by requiring 
the Gibbs conditions on the global conservation of the electric and the baryon charges. The relevance of small 
deviations from the standard extensive statistics is investigated in the context of intermediate and high energy 
heavy-ion collisions. 
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1. INTRODUCTION 

The determination of the properties of nuclear matter as 
functions of densities, temperature and protons fraction is a 
fundamental task in nuclear and subnuclear physics. Several 
experimental observations and theoretical calculations clearly 
indicate that hadrons dissociate into a plasma of their elemen- 
tary constituents, quarks and gluons (QGP), at density sev- 
eral times the nuclear matter density and/or at temperature 
above few hundreds MeV, which is the critical temperature 
T c of the transition from the QGP phase to the hadronic gas 
phase and viceversa. Such a QGP is expected to have oc- 
curred in the early stages of the Universe and can be found 
in dense and hot stars, neutron stars, nucleus-nucleus high en- 
ergy collisions where heavy ions are accelerated to relativistic 
energies [1]. After collision, a fireball is created which may 
realize the conditions of the QGP. The plasma then expands, 
cools, freezes-out into hadrons, photons, leptons that are de- 
tected and analyzed J2j] . 

It is a rather common opinion that, because of the ex- 
treme conditions of density and temperature in ultrarelativis- 
tic heavy ion collisions, memory effects and long-range color 
interactions give rise to the presence of non-Markovian pro- 
cesses in the kinetic equation affecting the thermalization pro- 
cess toward equilibrium as well as the standard equilibrium 
distribution JJ 0, Ht] . A rigorous determination of the con- 
ditions that produce a nonextensive behavior, due to memory 
effects and/or long-range interactions, should be based on mi- 
croscopic calculations relative to the parton plasma originated 
during the high energy collisions. At this stage we limit our- 
selves to consider the problem from a qualitative point of view 
on the basis of the existing theoretical calculations and exper- 
imental evidences. 

On the other hand, over the last years, there has been an 
increasing evidence that the generalized non-extensive statis- 
tical mechanics, proposed by Tsallis [g |8[1 and character- 
ized by a power-law stationary particle distribution, can be 
considered as a basis for a theoretical framework appropriate 
to incorporate, at least to some extent and without going into 
microscopic dynamical description, long-range interactions, 
long-range microscopic memories and/or fractal space-time 
constraints. A considerable variety of physical issues show a 
quantitative agreement between experimental data and theo- 
retical analysis based on Tsallis' thermostatistics. In partic- 



ular, there is a growing interest in high energy physics ap- 
plications of non-extensive statistics H [T(| Qjl El 0, Q. 
Several authors outline the possibility that experimental ob- 
servations in relativistic heavy-ion collisions can reflect non- 
extensive statistical mechanics effects during the early stage 
of the collisions and the thermalization evolution of the sys- 
tem H [H d [LI H [HIU • In this context, it is relevant 
to observe that the statistical origin of the nonextensive statis- 
tics lies in the deformation of the Boltzmann entropy. Let 
us note that in literature other statistical generalizations are 
present, such as the ^-deformed thermostatistics inspired to 
^-deformed quan tu m alg ebra and quantum groups with a dif- 
ferent origin im m sim . 

From the above considerations, it appears reasonable that 
in regime of high density and temperature both hadron and 
quark-gluon Equation Of State (EOS) can be sensibly affected 
by nonextensive statistical effects B25I1 . Furthermore, in this 
context it is very remarkable to observe that the relevance of 
these effects on the relativistic hadronic equation of state has 
also been recently investigated in Ref. lErill . 

The aim of this paper is to study the behavior of the nuclear 
equation of state at finite temperature and baryon density and 
to explore the existence of a hadron-quark mixed phase at a 
fixed value of the proton fraction Z/A. 



2. RELATIVISTIC NONEXTENSIVE 
THERMODYNAMICS 

In order to study, from a phenomenological point of view, 
experimental observable in relativistic heavy-ion collisions, in 
this Section we present the basic macroscopic thermodynamic 
variables and kinetic theory in the language of the nonexten- 
sive relativistic kinetic theory. 

Let us start by introducing the particle four-flow in the 
phase space as 12711 



A» J P 



and the energy-momentum flow as 



&q J p 



(1) 



(2) 
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where we have set h 



1, X = = (f,x), p = p p 



(p ,p), p° = vP^+ff? being the relativistic energy and 
f(x,p) the particle distribution function. In the above Z q = 
J d£l[f(x 1 p)] q is the nonextensive partition function, dQ. 
stands for the corresponding phase space volume element and 
q is the deformation parameter. The four- vector iV = («,j) 
represents the probability density n = n(x) (which is normal- 
ized to unity) and the probability flow j = j(jc). The energy- 
momentum tensor contains the normalized g-mean expecta- 
tion value of the energy density, as well as the energy flow, 
the momentum and the momentum flow per particle. Its ex- 
pression follows directly from the definition of the mean q- 
expectation value in nonextensive statistics ||7|]; for this reason 
7^ v it is given in terms of f q (x,p). 

On the basis of the above definitions, one can show that it is 
possible to obtain a generalized non-linear relativistic Boltz- 
mann equation 027I1 



in M \f(x,p)]* = C q (x,p) 



(3) 



where the function C q (x,p) implicitly defines a generalized 
nonextensive collision term 



1 [d 3 Pl d 3 p' d 3 p' l f . . , , 

Cq{x,p) = - / — jg JS-ikqlf J l]W(p ,pi\p,pi) 

LJ Pi p Pi 

-h q {f,fi}W(p, Pl \p',p\)}. (4) 



At equilibrium, the solution of the above Boltzmann equa- 
tion is a relativistic Tsallis-like (power law) distribution and 
can be written as 



feqip) = TT 



l-(l-q) 



P"Uu 



1/(1-9) 



(7) 



where is the hydrodynamic four- velocity [28] and f eq de- 
pends only on the momentum in the absence of an external 
field. At this stage, T is a free parameter and only in the 
derivation of the equation of state it will be identified with 
the physical temperature. 

We are able now to evaluate explicitly all other thermo- 
dynamic variables and provide a complete macroscopic de- 
scription of a relativistic system at the equilibrium. Con- 
sidering the decomposition of the energy-momentum tensor: 
P ,v = sU^U* —PA^, where e is the energy density, P the 
pressure and = g^ v — U P U V , the equilibrium pressure can 
be calculated as 



1 



1 



d 3 p 



P = --T^A, V = -— / -i-pVA^fi^) • (8) 



Setting T = p /T and z~m/T, the above integral can be ex- 
pressed as 



4n 



m 2 T 2 K 2 (q,z) 



(9) 



Here W(p 1 p\ \p',p'i) is the transition rate between a two- 
particle state with initial four-momenta p and p\ and a fi- 
nal state with four-momenta p' and p\\ h q [f,f\\ is the q- 
correlation function relative to two particles in the same space- 
time position but with different four- momenta p and pi, re- 
spectively. Such a transport equation conserves the probabil- 
ity normalization (number of particles) and is consistent with 
the energy-momentum conservation laws in the framework of 
the normalized g-mean expectation value. Moreover, the col- 
lision term contains a generalized expression of the molecular 
chaos and for q > implies the validity of a generalized H- 
theorem, if the following, nonextensive, local four-density en- 
tropy is assumed (henceforward we shall set Boltzmann con- 
stant k B to unity) 



d 3 P 



pV[(*,/0]*H/(*>/0-i] 



(5) 



the inverse function of the Tsallis' ^-exponential function 



where we have used the definition ln^x 



e q (x) = [l + {l-q)x] 



i/(i-?) 



(6) 



which satisfies the property e q (ln q x) = x. 

The above expression is written in a covariant form, in fact 
S q = (S q ,S' q ), with i = 1,2,3, and correctly transforms as a 
four- vector under Lorentz transformations B27I1 . where S q is 
the relativistic Tsallis nonextensive local entropy density and 
S'q is the Tsallis entropy flow. Note that for q —> 1, Eqj 
reduces to the well known four-flow entropy expression 



where we have introduced the ^-modified Bessel function of 
the second kind as follows 



K„ 



2"n\ 1 



dx(x 2 -z 2 y 



(10) 



and e q (x) is the ^-modified exponential defined in Eq.©. 

Similarly, the energy density e can be obtained from the 
following expression 



e = T^U M U y = 



1 



d 3 p 



(p"U,) 2 fU P ), (11) 



and, after performing the integration, it can be cast into the 
compact expression: 



4n 



K 2 (q,z) , Ki(q,z) 



Thus the energy per particle e — e/n is 

Ki(q,z) 



' — 3 T + m 



K 2 (q,z) ' 



(12) 



(13) 



which has the same structure of the relativistic expression ob- 
tained in the framework of the equilibrium Boltzmann-Gibbs 
statistics JH]. 

In the non-relativistic limit (pel) the energy per particle 
reduces to the well-known expression 



(14) 
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and no explicit ^-dependence is left over. 

Hence from the above results it appears that, in search- 
ing for the relevance of nonextensive statistical effects, both 
microscopic observable, such as particle distribution, correla- 
tion functions, fluctuations of thermodynamical variables, and 
macroscopic variables, such as energy density or pressure, can 
be affected by the deformation parameter q. 

In this context, it appears relevant to observe that, in 
Ref. Hill nonextensive Boltzmann equation has been studied 
and proposed for describing the hadronization of quark mat- 
ter. Moreover, starting from the above generalized relativistic 
kinetic equations, in Ref. [18] the authors have recently for- 
mulated a nonextensive hydrodynamic model for multiparti- 
cle production processes in relativistic heavy-ion collisions. 
These works represent an important bridge for a close con- 
nection between a microscopic nonextensive model and ex- 
perimental observable. 

Finally, let us remind the reader that for a system of par- 
ticles in a degenerate regime the above classical distribution 
function © has to be modified by including the fermion and 
boson quantum statistical prescriptions. For a dilute gas of 
particles and/or for small deviations from the standard exten- 
sive statistics (q ss 1) the equilibrium distribution function, in 
the grand canonical ensemble, can be written as 1 29] 

" (M = [i + ( q - 1 )(£(*) -AO/r] 1 /^-!) ± i ' (15) 

where the sign + stands for fermions and — for bosons: hence 
all previous results can be easily extended to the case of quan- 
tum statistical mechanics. 



collisions is that, for symmetric or nearly symmetric nuclear 
matter, the critical density (at low temperatures) appears to be 
considerably larger than nuclear matter saturation density po- 
Concerning non-symmetric matter, general arguments based 
on Pauli principle suggest that the critical density decreases 
with Z/A. Therefore, the transition's critical densities are ex- 
pected to sensibly depend on the isospin of the system lO . 
Moreover, the analysis of observations of neutron stars, which 
are composed of P-stable matter for which Z/A < 0.1 (the 
matter constituting neutron stars is strongly isospin asymmet- 
ric, being composed of a large amount of neutrons and a small 
fraction of protons) can also provide hints on the structure of 
extremely asymmetric matter at high density. No data on the 
quark deconfinement transition are at the moment available 
for intermediate values of Z/A. Recently, it has been proposed 
by several groups to produce unstable neutron-rich beams at 
intermediate energies. These new experiments can open the 
possibility to explore in laboratory the isospin dependence of 
the critical densities. 

The scenario we are going to explore in this last Section 
corresponds to the situation realized in experiments at not 
too high energy. In this condition, only a small fraction of 
strangeness can be produced and, therefore, we limit our- 
selves to study the deconfinement transition from nucleonic 
matter into up and down quark matter. In the next two subsec- 
tions, we will study the two corresponding EOSs separately, 
on the basis of the previously reported nonextensive relativis- 
tic thermodynamic relations. The existence of the hadron- 
quark mixed phase will be studied in the third subsection. This 
investigation may be helpful also in view of the future experi- 
ments planned, e.g., at the facility FAIR at GSI Q. 



3. NONEXTENSIVE NUCLEAR EQUATION OF STATE 

The EOS at densities below the saturation density of nu- 
clear matter po « 0.16 fm~ 3 is relatively well known due to 
the large amount of experimental nuclear data available. At 
larger densities there are many uncertainties due the lack of 
experimental data; the strong repulsion at short distances of 
nuclear force makes, in fact, the compression of nuclear mat- 
ter to larger densities quite difficult. Otherwise, in relativistic 
heavy ion collision the baryon densities can reach values of a 
few times po and the temperature can exceed the deconfine- 
ment critical temperature T c w 170 MeV. 

As partially discussed in the Introduction, hadronic matter 
is expected to undergo a phase transition into a deconfined 
phase of quarks and gluons at large densities and/or high tem- 
peratures. However, the extraction of experimental informa- 
tion about the EOS of matter at large densities and tempera- 
tures at intermediate and high energy heavy-ion collisions is 
very complicated. Possible indirect indications of a softening 
of the EOS at the energies reached at AGS have been dis- 
cussed several times in the literature iTsol l3lll . In particular, 
a recent analysis [32] based on a 3-fluid dynamics simulation 
suggests a progressive softening of the EOS tested through 
heavy-ion collisions at energies ranging from 2A GeV up to 
8A GeV. On the other hand, the information coming from ex- 
periments with heavy-ions at intermediate and high energy 



3.1. Nonextensive hadronic equation of state 

The relativistic, field theoretical approach to nuclear EOS 
was used first by Walecka and Boguta-Bodmer in the mid- 
1970s 11351, 13611 . This theory describes the interaction between 
nucleons through the exchange of two mesons, the scalar field 
a and the vector field co. The model of Walecka has two free 
parameters: the two ratios between the nucleon-meson cou- 
pling constants and the masses of the mesons. The saturation 
density and binding energy per nucleon (calculated at the the 
saturation density) of nuclear matter can be fitted exactly in 
the simplest version of this model but other properties of nu- 
clear matter, as e.g. incompressibility, cannot be reproduced. 
To overcome these difficulties, the model has been modified 
introducing in the Lagrangian two terms of self-interaction for 
the a which are crucial to reproduce the empirical incompress- 
ibility of nuclear matter and the effective mass of nucleons M* 
(again calculated at the saturation density). Moreover, the in- 
troduction of an isovector meson p allows to reproduce the 
correct value of the empirical symmetry energy 11371 13811 . 

In the following, we will use a relativistic mean field self- 
consistent theory of nuclear matter in which nucleons interact 
through the nuclear force mediated by the exchange of virtual 
isoscalar and isovector mesons (a,co,p) with a Lagrangian 
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density [3 



\p^--F^--G^&\, 



(16) 



where M — 939 MeV is the vacuum baryon mass, (0,tt>u) 
are the isoscalar (scalar,vector) meson fields, while p^ is the 
corresponding iso vector ones. The field strength tensors for 
the vector mesons are given by the usual expressions = 
d/j(i3 v — d v (x) p , G pV = d p p v — dvP/j' an d the U(g) is a nonlinear 
potential of o meson 



U(a) = iflo 3 



1 



-bo* 



(17) 



This last term is usually introduced to achieve a reasonable 
compression modulus for equilibrium nuclear matter. 
The field equations in a mean field approximation are 

(iy^ - (Af - g a o) - g a r°^o - gpYVpcOv = , 



m\(S- 



-ao 2 + bo i =g a < W>=gaPs, 



m^COo = gm < vYV >= gapB , 
mjpo =gp< VfnfViV >= SpPi , 



(18) 



where pj = p p — p„ is the isospin density, and p$ are the 
baryon and the scalar densities, respectively. They are given 
by 



i=n.p 



i— «,p 



(2jt) 
d 3 k M 



(2n 



3 +^(*)]. 



(19) 



where rij(k) and n, (£) are the g'-deformed fermion particle and 
antiparticle distributions given in Eq.(fT3]>; more explicitly 



ni{k) = 



m(k) = 



[l + (q-l)(Ef(k)-rf)/T]V<.*-V + l : 
1 

[1 + (q - l){E*(k) + jU *)/r] !/(<?-!) + 1 



(20) 



(21) 



The nucleon effective energy is defined as E*(k) = 
\/k 2 +M* 2 , where M* = M, - g a a. The effective chemi- 
cal potentials /j* are given in terms of the vector meson mean 
fields fj* = /jj — g m 03o T^pPo (— proton, + neutron), where 
fjj are the thermodynamical chemical potentials /j, = 3e/3p,. 
At zero temperature they reduce to the Fermi energies Ep ; = 

kj^f and the nonextensive statistical effects disappear. 
The meson fields (a, COo and po) are obtained as a solution of 
the field equations in mean field approximation and the related 
couplings (g a , g m and g p ) are the free parameter of the model 
& 



P [Gev/fm J ] 
0.4 r 



q=1.2_ 



0.9 



1 1.1 1.2 1.3 



/j[GeV] 



FIG. 1: Hadronic equation of state: pressure versus baryon chemical 
potential for different values of q. In the figure T = 100 MeV and 
Z/A = 0.4. 



On the basis of Eqs.(f2|l, ([8]) and ( fTTl i. the pressure and the 
energy density can be written as 



^m 2 a a 2 - U(a) + ijf^cog + \m\p\ , 



e = 



2E 



d 3 k 



(27T)3 



3 £f(fc)[nf(fe)+nf(fc)] 
»4o 2 + U(a) + l -ml&l + l -m\p\ , 



(22) 



(23) 



Note that statistical mechanics enters as an external ingre- 
dient in the functional form of the "free" particle distribution 
of Eq. $15[ . Since all the equations must be solved in a self- 
consistent way, the presence of nonextensive statistical effects 
in the particle distribution function influences the many-body 
interaction in the mean field self-consistent solutions obtained 
for the meson fields. 

In Fig. [T| we report the resulting hadronic EOS: pressure as 
a function of baryon chemical potential for different values of 
q. Since in the previous investigations we have phenomeno- 
logically obtained values of q greater than unity |5l |20[] , we 
will concentrate our analysis to q > 1 . The results are plotted 
at the temperature T = 100 MeV, at fixed value of Z/A = 0.4 
and we have used the GM2 set of parameters of Ref. lBal . The 
range of the considered baryon density and the chosen values 
of the parameters correspond to a physical situation which can 
be realized in the recently proposed high energy heavy-ion 
collisions experiment at GSI IB9I1 . 



3.2. Nonextensive QGP equation of state 

In high density hadronic matter, baryons are forced to stay 
so close one to another that they would overlap. At such densi- 
ties, constituent quarks are shared by neighboring baryons and 
should eventually become mobile over a distance larger than 
the typical size of one baryon. This means that quarks become 
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deconfined and that at large densities and/or high tempera- 
tures they are the real degrees of freedom of strongly interact- 
ing matter instead of baryons. The process of deconfinement 
and the EOS of quark matter can in principle be described by 
quantum chromodynamics. However, in energy density range 
reached in the core of compact stars and in relativistic heavy- 
ion collisions, non-perturbative effects in the complex theory 
of QCD are present. For these reasons, simple phenomeno- 
logical models are usually adopted to describe quark matter 
ESE110. 

In this study, due to its simplicity, we use the MIT bag 
model [40]. In this model, quark matter is described as a gas 
of free quarks with massless up and down quarks. All the non- 
perturbative effects are simulated by the bag constant B which 
represents the pressure of the vacuum. Following this line, 
the pressure, energy density and baryon number density for a 
relativistic Fermi gas of quarks in the framework of nonexten- 
sive statistics (see Eqs.(Q]i, ©, © and ( fTTT i) can be written, 
respectively, as 



iVk k l[n<* f {k) + rf f (k)]-B, (24) 



If I r^e f [n}(k)+jf f (k)]+B, (25) 
/=«,</ ^ ' 



If 



d 3 k 



3 /~> (27!) 



[n f (k)-n f (k)], 



(26) 



where the quark degeneracy for each flavor is jf = 6, e/ = 
{k 2 +m 2 f y/ 2 ,n f (k) <mdn f (k) are the ^-deformed particle and 
antiparticle quark distributions 



P [Gev/fm J ] 



q=1.2 



:.2 



1 1.1 1.2 1.3 



FIG. 2: The same as in Fig.fT]for the case of the quark-gluon equation 
of state. 



As previously discussed, in this investigation we are lim- 
iting our study to the two-flavor (/ = u, d) massless quarks. 
As already remarked, this appears rather well justified for the 
application to heavy ion collisions at relativistic (but not ultra- 
relativistic) energies, the fraction of strangeness produced at 
these energies being small l43l 14411 . 

In Fig. |2j we report the pressure as a function of the baryon 
chemical potential for massless quarks u, d and gluons, for 
different values of q. As in Fig.[T] the results are plotted at the 
temperature T = 100 MeV and at a fixed value of Z/A = 0.4; 
the bag parameter is B 1 / 4 =170 MeV. In both figures Q] and [2] 
one can observe sizable effects in the behavior of the EOS 
even for small deviations from the standard statistics. 



n f (k) 



n f (k) 



[i + (q- 1 )(>/(£) -^/Vr] 1 /^/- 1 ) + 1 ' 

1 

[1 + (q - l)(e f (k) + W )/r] 1 /(?-l) + 1 ' 



(27) 



(28) 



Similar expressions for the pressure and the energy density 
can be written for the gluons treating them as a massless q- 
deformed Bose gas with zero chemical potential. Explicitly, 
we can calculate the nonextensive pressure P g and density en- 
ergy Eg for gluons as 



J g f 00 d 3 k k 
8 ~JJo (27t) 3 [\ + {q-\)k/T]<il(l- l )-\ 
E g =3P g , 



, (29) 
(30) 



with the gluon degeneracy factor y g = 16. Let us observe that, 
in the limit q — » 1, we recover the usual analytical expression 
for the pressure of gluons: P g = 87t 2 /45 T 4 . 

Since one has to employ the fermion (boson) nonextensive 
distributions, the results are not analytical, even in the mass- 
less quark approximation. Hence a numerical evaluations of 
the integrals in Eq.s (I24li-(|26]| and Eq.s (f29b — (f30b must be 
performed. Let us remember that a similar calculation, only 
for the quark-gluon phase, was also performed in Ref. 14511 by 
studying the phase transition diagram. 



3.3. Mixed hadron-quark phase 

In this subsection we investigate the hadron-quark phase 
transition at finite temperature and baryon chemical potential 
by means of the previous relativistic EOSs. Lattice calcula- 
tions predict a critical phase transition temperature T c of about 
170 MeV, corresponding to a critical energy density e c « 1 
GeV/fm 3 (]]]. In a theory with only gluons and no quarks, the 
transition turns out to be of first order. In nature, since the u 
and d quarks have a small mass, while the strange quark has a 
somewhat larger mass, the phase transition is predicted to be 
a smooth cross over. However, since it occurs over a very nar- 
row range of temperatures, the transition, for several practical 
purposes, can still be considered of first order. Indeed the lat- 
tice data with 2 or 3 dynamical flavours are not precise enough 
to unambigously control the difference between the two situ- 
ations. Thus, by considering the deconfinement transition at 
finite density as a the first order one, a mixed phase can be 
formed, which is typically described using the two separate 
equations of state, one for the hadronic and one for the quark 
phase. 

To describe the mixed phase we use the Gibbs formalism, 
which in Ref. i46ll has been applied to systems where more 
than one conserved charge is present. In this contribution we 
are studying the formation of a mixed phase in which both 
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baryon number and isospin charge are preserved. The main 
result of this formalism is that, at variance with the so-called 
Maxwell construction, the pressure in the mixed phase is not 
constant and therefore the nuclear incompressibility does not 
vanish. It is important to notice that from the viewpoint of 
Ehrenfest's definition, a phase transition with two conserved 
charges is considered, in the literature, not of first, but of sec- 
ond order H3] . 

The structure of the mixed phase is obtained by imposing 
the Gibbs conditions for chemical potentials and pressure and 
by requiring the global conservation of the total baryon (B) 
and isospin densities (I) in the hadronic phase (H) and in the 
quark phase (Q) 

(H) (Q) 
(H) (Q) 



P« = (l-X)pf + XPb, 
P/ = (1-X)pf + XP?, 



P [Gev/fm J ] 

0.25 

0.2 
0.15 

0.1 
0.05 



q=i-i.-- 



1.05 1.1 1.15 1.2 1.25 1.3 1.35 



fJ[GeV] 



FIG. 3: Pressure versus baryon chemical potential in the mixed 
hadron-quark phase for different values of q. The long-dashed line 
corresponds to q = 1.05. The temperature is fixed a T=60 MeV and 
the proton fraction Z/A = 0.4. 



(31) 



where % is the fraction of quark matter in the mixed phase. 
In this way we can obtain the binodal surface which gives the 
phase coexistence region in the (T, Pb,P/) space. For a fixed 
value of the conserved charge p/, related to the proton fraction 
Z/A = (l+p//pfl)/2, we study the boundaries of the mixed 
phase region in the (T,ps) plane. We are particularly inter- 
ested in the lower baryon density (or baryon chemical poten- 
tial) border, i.e. the critical/transition density p cr , in order to 
check the possibility of reaching such (r,p CT ,p/) conditions 
in a transient state during a heavy-ion collision at relativistic 
energies. 

In Fig. [3] we report the pressure versus the baryon chemical 
potential, in Fig. [4] the pressure as a function of the energy 
density and in Fig. [5] we report the pressure versus baryon 
density (in units of the nuclear saturation density po) in the 
mixed hadron-quark phase for different values of q. For the 
hadronic phase we have used the so-called GM2 set of param- 
eters l38ll and in the quark phase the bag parameter is fixed 
to B 1 I A =\1Q MeV. The temperature is fixed at T = 60 MeV 
and the proton fraction at Z/A=0.4, physical values which are 
estimated to be realistic for high energy heavy-ion collisions. 
The mixed hadron-quark phase starts at p = 3.75 po for q=l, 
atp = 3.31 po for q = 1.05 and at p = 2.72po for q = 1.1. It is 
important to observe that for q = 1.1 the second critical transi- 
tion density is also reached, separating the mixed phase from 
the pure quark-gluon matter phase, at p = 4.29 po while for 
q = 1.05 the second critical density is reached at p = 5.0 po 
and at p = 5.57po for q=l. 

As a concluding remark we note that nonextensive statisti- 
cal effects become extremely relevant at large baryon density 
and energy density, as the ones which can be reached in high 
energy collisions experiments. This fact can be an important 
ingredient in the realization of a hydrodynamic model as well 
as to obtain a deeper microscopic connection with the experi- 
mental observable. 
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FIG. 4: Pressure versus energy density in the mixed hadron-quark 
phase for different values of q. The temperature is fixed a T=60 MeV 
and the proton fraction Z/A = 0.4. 
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FIG. 5: Pressure versus baryon density (in units of the nuclear satura- 
tion density po) in the mixed hadron-quark phase for different values 
of q. 
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